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It is shown that, while Wigner and Liouville functions transform in an identical way under 
linear symplectic maps, in general they do not transform identically for nonlinear symplectic 
maps. Instead there are "quantum corrections" whose ?i ^ limit may be very complicated. 
Examples of the behavior of Wigner functions in the h limit are given in order to 
examine to what extent the corresponding Liouville densities are recovered. 
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I. INTRODUCTION 



Consider an ensemble of non- interacting particles. 
Classically such an ensemble may be described by a Li- 
ouville density in phase space. Denote the collection of 
phase-space variables by the symbol z, 



z = (gi,g2,93;Pi,P2,P3)- 



(1.1) 



Let w{z) be the Liouville density at the point z. Then, 
by definition, the number of particles dPN in the phase- 
space volume dPz is given by the relation 



d^N = w{z)(fz. 



(1.2) 



Suppose the ensemble of particles is initially character- 
ized by a distribution function . Let be a symplectic 
map (canonical transformation) that sends initial points 

to final points z^ , 



(1.3) 



This map has the associated Jacobian matrix M{z^) de- 
fined by the relation 



Mal,(z') = dzi/dz 



(1.4) 



By the definition of M being a symplectic map, the ma- 
trix M must be a symplectic matrix. 



M^JM = J. 



(1.5) 



Here J is the usual fundamental matrix having the 3x3 
block form 

It can be shown from (1.5) that M has unit determinant 



detM = 1. 



(1.7) 



Correspondingly, if is the final distribution function 
resulting from the action of AA, Liouville's theorem fol- 
lows as a consequence of (1.7) and (1.2), 



wf{zf)^w\z'^). (1.8) 

Now insert (1.3) into (1.8) to find the relation 

w^{zf) = w\M-^z^). (1.9) 

Since z^ is a generic point, the relation (1.9) may equally 
well be written in the form 



wf{z) = w\M-^z). 



(1.10) 



This equation describes how the Liouville function trans- 
forms under the action of a general symplectic map. 

Next consider a quantum description of this same en- 
semble of noninteracting particles. In the quantum de- 
scription, an ensemble is characterized by a density oper- 
ator p. Associated with any density operator is a Wigner 
function M^(z) that in some ways is the quantum analog 
of the Liouville function w{z) j^. Also, there are unitary 
transformations lA that are the quantum analog of the 
symplectic maps M. It is therefore tempting to consider 
the possibility 



Wf{z) = W'iM-'^z). 



(1.11) 



It is known that (1.11) is true in the case that is a 
translation or a linear transformation The purpose 
of this note is to show, by a simple counter-example, 
that (1.11) is not in general true for a nonlinear sym- 
plectic map. Section 1.2 describes needed Lie-algebraic 
tools in the classical case. Section 1.3 describes quantum 
concepts and defines the Wigner function. Section 1.4 
presents the counter-example, and discusses the h ^ 
limit for sample Wigner functions. A final section dis- 
cusses various aspects of what is already known and avail- 
able in the literature about Wigner functions, and shows 
that in this context the failure of (1.11) in the nonlinear 
case is no surprise. 

Why should transformation properties of the Wigner 
function be of interest? 

First, there is the general question of the relation be- 
tween quantum and classical mechanics and the way that 
quantum mechanics becomes classical in the limit ?i — > 0. 
Second, a light beam may be described, using ray theory. 
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in terms of a Liouville function. In wave theory, where 
the "reduced" wavelength ^ plays a role analogous to 
Planck's constant, the same beam may be described by 
a Wigner function. In the ray description, aberration ef- 
fects are characterized by nonlinear symplectic maps . 
Correspondingly, in a wave description, one would like to 
know the effect of aberrations on the Wigner function. 

II. LIE ALGEBRAIC TOOLS 

Let /(z) be any function of the phase-space variables 
z. Associated with / is a dijferential operator, called a 
Lie operator, that will be denoted by the symbol :/: and 
is defined by the rule 

■■f--= J2(df/dq.)id/dp.) - (2.1) 

i 

When :/: acts on any other phase-space function g(z), 
the result is a Poisson bracket, 

:/:.?=[/, 5]- (2.2) 

Powers of Lie operators are defined by repeated appli- 
cation, with :/:" defined to be the identity operator, 

:/:° 9 = 9, 
■■f-' 9^[f,9] 

■.f:'g = :f: (:/: .9) = [/,[/, .g]], etc. (2.3) 

Once powers have been defined, power series are also de- 
fined. Of particular interest is the exponential series, 

00 

exp(:/:) = e-f- = ^ :/:^ (2.4) 

This series is called a Lie transformation. According to 
(2.3) and (2.4), a Lie transformation acts on the function 
g{z) to give the result 

exp(:/:)g = g + [/, ff] + [/, [/, <?]]/2! + • • • . (2.5) 

Suppose / is any function of the initial phase-space 
variables z*. Define final variables z-^ by the rule 

zf - cxp(:/(z'):)z^ (2.6) 

Comparison of (1.3) and (2.6) indicates that in this case 
the map A4 can be written in the form 

M = exp(:/:). (2.7) 

It can be shown that any map of the form (2.7) is a 
symplectic map. Conversely, any symplectic map can be 
written as a product of Lie transformations jl] . 

Three simple Lie transformations are of particular in- 
terest for subsequent use. All apply to the case of a 
2-dimensional phase space described by the variables q 



and p. First consider the symplectic map given by (2.7) 
when / has the form 

f{z) = aq. (2.8) 

Here a is a parameter. Then one finds from (2.3) and 
(2.4) the results 

q^ ^ Mq'' = exp[:ag*:]g* = q\ 

pi = Mp^ = exp[:Q;g*:]p* = p' + a. (2.9) 

Second, suppose / is of the form 

/(z) - (2.10) 

Then one finds the results 

qf = Mq' = cxp[:a(g')V2:]g' = q\ 

pf = Mp' = exp[:a{q'f /2:]p' = p' + aq\ (2.11) 

Finally, suppose / is of the form 

/(z) = aq^/i. (2.12) 

Then one finds the results 

qS ^ Mq' = exp[:a(g^)V3:]<7' = q\ 

pf ^ Mp' ^ exp[:a{q' f/3:]p' ^ p' + a(q' f. (2.13) 

Transformation (2.9) is a translation, (2.11) is a linear 
transformation, and (2.13) is a nonlinear transformation. 

III. QUANTUM CONCEPTS 

For simplicity of presentation, the quantum treatment 
will be described for the case of a 2-dimensional phase 
space and the associated operators Q and P. General- 
ization to higher dimensions will be obvious. 

Suppose A is any operator (presumably involving Q 
and P). Associated with this operator is the function 
yl(g,p), called the Weyl transform of A, and defined by 
the rule 

A{q,p)^ J dy{q + y/2\A\q~y/2)expi-ipy/h). (3.1) 

Here \q), as usual, denotes a position eigenstate of the 
operator Q. Conversely, the underlying operator A can 
be recovered from a knowledge of A by using an inverse 
Weyl transformation, 

A{Q,P) = {2Trn)-^ J d(t)dTdqdpA{q,p) 

X exp{(V?i)[0(Q -q)+ t{P - p)]}. (3.2) 
The reader should verify that if 

A^T{P) + V{Q), (3.3) 
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then, by the Weyl correspondence, 

A = T{p) + V{q), (3.4) 

and vice versa. FinaUy, the Weyl transform has the prop- 
erty that, if A and B are any two operators, there is the 
relation 

{2T:h)-^ j dqdpA{q,p)B{q,p) = tr{AB). (3.5) 

Consider a quantum ensemble consisting of N non- 
interacting members. Let the j'th member of the en- 
semble be described by the state vector \tpj). Then the 
density operator p for this ensemble is defined by the rule 

N 

p^{l/N)J2\^,){i^j\- (3.6) 
i=i 

Evidently by construction the density operator is Her- 
mitian and positive definite, and has unit trace. If A is 
any operator, its expectation value for the j'th member 
in the ensemble is (tpjlAl^j) , and the ensemble average 
of this expectation value is 

N 

mAmu = (i/iv)^(v.,iii^,). (3.7) 

Let \xn) be any orthonormal set of basis vectors. The 
respresentation of the identity / in this basis is given by 
the relation 

n 

Now insert (3.8) into (3.7) to find the familiar formula 
that relates a combined quantum-mechanical and ensem- 
ble average to a trace, 

N 

[{^pM\^,)U = ii/N)j2{i^,\A\^,) 

N 

= (l/iV)^(V;,|i/|^,) 

N 

^{l/N)Y,J2(^j\A\Xn){XnHj) 
n j = l 

= J2{Xn\{l/N)Y,\^,){^^,\A\Xn) 

n j 
n 

The Wigner function associated with any density oper- 
ator p is defined to be the Weyl transform of the operator 
p/i2nn), 

Wiq,p) = {2^h)-^ j dy{q + y/2\p\q - y/2) exp{-ipy/h). 

(3.10) 



It follows from this definition and (3.5) that W has the 
properties 

J dqdpWiq,p)Aiq,p) = tr{pA), (3.11) 

j dqdpW{q,p) = tr{pl) = tr{p) = 1. (3.12) 

The left side of (3.11) resembles a classical phase-space 
average, and the right side of (3.11) is a quantum- 
mechanical ensemble average. Thus, in this respect, 
the Wigner function W resembles a normalized Liouville 
function. It can also be verified that 

j dpW{q,p) = {q\p\q) > 0, (3.13) 

j dqW{q,p)^{p\p\p)>0, (3.14) 

SO that any partially integrated Wigner function behaves 
like a position or momentum density. However, unlike 
a Liouville density, in general the Wigner function itself 
may sometimes be negative |^]. 

IV. SIMPLE COUNTER-EXAMPLE 

Suppose that each member of a quantum ensemble 
is acted upon by a common unitary transformation U. 
Then the state vector for the j'th member of the 
ensemble becomes \'ip'j) with 

l^;)=wiv,), (4.1) 

and correspondingly, according to (3.6), the density op- 
erator p for the ensemble becomes p' with 

p'^UpU-^. (4.2) 

It follows from the definition (3.10) that the transformed 
Wigner function W' associated with p' is given by the 
relation 

W\q,p) = (27r?i)-l j dy{q + y/2\Up U-^\q ~ y/2) 

X exp(-ipy/?i). (4.3) 

In analogy with (2.7) and (2.8), suppose U is of the 
form 

U = exp(iaQ/;i). (4.4) 

(This analogy is motivated by the close resemblance be- 
tween the Poisson bracket Lie algebra of Classical Me- 
chanics and the commutator Lie algebra of Quantum Me- 
chanics.) For a position eigenstate \q) there is the result 
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U\q) = exp(iQ;Q/?i)|(7) = ex-]^{iaq/%)\q) , (4.5) 
and for the momentum eigenstate \p) there is the result 
U\p) = exp(mQ/?i)|p) = \p + a). (4.6) 
This latter result follows from the fundamental relation 
\p) = (27r?i)-i/2 f dqeM^q/m), (4.7) 



or, equivalently, from the fundamental commutation rule 
{Q,P}^QP-PQ^ihI. (4.8) 



U-^lq - y/2) = cxp[-taiq - v/2f/{2n)]\q - y/2), 

(4.16) 

and the "exponent addition" relation analogous to (4.12) 
takes the form 

[-ia{q - y/2f/{2n)] + [+ia{q + y/2f/{2h)] + {-ipy/K) 
= -i{p-aq)y/h. (4.17) 

The reader is urged to verify this result in the solitude of 
some evening. It follows that in this case 



Evidently, (4.5) and (4.6) are the quantum analog of 
(2.9). 

What is W when U is given by (4.4)? Examine the in- 
gredients of the integrand for (4.3). From (4.5) it follows 
that 



W\q,p) = W{q,p-aq), 

and consequently (1.11) again holds. 
Finally, suppose that U is of the form 



(4.18) 



W = exp[iaQV(3?i)]- (4.19) 
In this case the exponent addition relation is 
U~^\q - y/2) = exp[-ia{q - y/2)/h]\q - y/2), (4.9) [_^^(^ _ y/2)3/(3?i)] + [+ia(g + y/2)V(3?i)] + {-tpy/h) 

^ -i{p-aq^)y/h + iay^/{l2h). (4.20) 

Correspondingly, W' is given by the relation 



{q + y/2\U - exp[+ia(g + y/2)/h]{q + y/2\. (4.10) 
Hence, (4.3) becomes 



W\q,p) = {2nn)-^ / dy{q + y/2\p\q-y/2) 



X exp[— ia(q — y/2)/h] exp[+ia(<? + y/2)/h] 
X exp{—ipy/h). (4-11) 

Moreover, addition of the exponents in (4.11) yields 

[-ia{q ~ y/2)/h] + [+ia{q + y/2)/h] + {-ipy/K) 
= -i{p-a)y/h. (4.12) 

Consequently, (4.3) takes the final form 

W'iq,p) = {27rh)-' f dy{q + y/2\p\q-y/2) 



(4.13) 



W'{q,p) = (2^n)-l J dy{q + y/2\p\q - y/2) 

X exp[~i(p - aq'^)y/h] ejip[iay^ / {12h)]. (4.21) 

Evidently (1.11) now fails to hold due to the offensive 
term exp[iay^ / {12'h)] in the integrand of (4.21). 

At this point we remark that a similar calcula- 
tion shows that (1.11) fails for any U of the form 
eicpliaQ" / (nh)] with n > 3. Moreover, it can be shown 
that any symplectic map can be written as a product of 
linear transformations and transformations of the form 
exp(a :g": /n) [|. If follows that (1.11) fails for all 
nonlinear transformations. 

What more can be said? Suppose the offensive term is 
expanded in a Taylor series, 



X exp[— i(p — a)y/h] 
= W{q,p-a). 

From (2.9) it follows that 

M^^q ~ exp(— ■.aq:)q — q, 

A4~^p = exp(— ■.aq:)p = p — a. 



exp[m2/3/(12?i)] = 1 -I- iay^/{l2n) 



(4.22) 



(4.14) 



Therefore, inspection of (4.13) and (4.14) illustrates that 
(1.11) holds for a translation. 
Next suppose that U is of the form 



U = exp[iaQV(2?i)] 



(4.15) 



The calculation proceeds as before. For example, there 
is the relation 



Insertion of this expansion into the integral (4.21) gives 
the result 

W'{q,p) = W{q,p-aq^) 

+ia/[{2^fi){l2h)] j dy{q + y/2\p\q~y/2) 

X exp[-i(p - aq^)y/h] [y^ + ■ ■ ■] . (4.23) 

Powers of y in the integrand can be obtained by differen- 
tiating the factor exjp{—ipy /h) with respect to p. Conse- 
quently, (4.23) can be recast in the form 

W'{q,p) 

= W{q,p - aq^) -|- aih^ / 12) {d / dp)^W {q,p - aq^) + ■■■. 

(4.24) 
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It is tempting to view (4.24) as a quantum correction to 
(1.11), and to speculate that (1-11) becomes exact in the 
limit h 0. This may be true in some cases providing 
W itself is well behaved as ?i ^ and remains sufficiently 
smooth so that the derivatives occurring in (4.24) remain 
bounded. For example, consider the simple harmonic 
oscillator described by the quantum Hamiltonian 

H = PV(2m) + (mwV2)g2 (4.25) 

corresponding to the classical Hamiltonian 

H{q,p) = p^/{2m) + {mw^/2)q^. (4.26) 

The density operator p for a canonical ensemble of such 
oscillators at temperature T, with /3 = l/{kT), is given 
by the relation 

p=[exp{-f3H)]/tr[exp{~(3H)]. (4.27) 

In this case the Wigner function is 

W{q,p) = [l/{7:h)]taiih{f3hLO /2)exp{-[2 /{huj)] 

xtanh{f3hLj/2)H{q,p)}. (4.28) 

Inspection of (4.28) shows that in this case W is well 
defined and smooth as ?i — > 0. Indeed, in this limit W 
becomes the classical Boltzmann distribution (Liouville 
function) , 

lim Wiq,p) = [/3c^/(2^)] eM~PH{q,p)]. (4.29) 

However, in many cases the Wigner function becomes 
ever more highly oscillatory as ?i ^ 0. Such behavior 
is typical of Wigner functions associated with ensembles 
that are relatively pure 0. For example, consider a per- 
fectly pure ensemble consisting of harmonic oscillators 
that are all in the same energy eigenstate. For such a 
pure ensemble, whose members are all in the eigenstate 
\n) for which the energy has the value 

{n+l/2)hcj, (4.30) 

the Wigner function is given by the equation Q 

W{n;q,p) = [{-1}" / {ttTi)] exp[-2H{q,p) / {Hu;)] 

xL^^)[4H{q,p)/{hcj)]. (4.31) 

Here lI""* is a Laguerre polynomial. 

For simplicity of presentation, suppose units are se- 
lected in such a way that both the spring constant k and 
the mass m have unit value. With this choice of units, 
the classical Hamiltonian becomes 

H={p' + q'-)/2, (4.32) 

and W takes the form 

W{n-r) = [{-ir/{TTn)]exp{-ryTi)L^°\2ryn). (4.33) 



Here we have introduced a radius r in phase space by the 
definition 

r^=p^ + q'^. (4.34) 

Inspection of (4.33) shows that W has an essential sin- 
gularity in 7i at ?i = 0. Moreover, differentiation oper- 
ators, such as (d/dp), produce powers of (l/K)- There- 
fore, since (4.24) contains (d/dp)^ and higher derivatives, 
what appear to be quantum corrections actually diverge 
'AS Ti ^ like {l/h) and higher powers of (1/^). Direct 
evaluation of (4.33) gives the result 

lim Ty(n;r) = if r > 0. (4.35) 

Moreover, as can be checked directly, (3.12) also holds. 
Therefore, in the limit, W{n; r) becomes a distribution 
that behaves like the delta function distribution centered 
on the phase-space origin, 

lim W{n, r) = [l/{2nr)]S{r) = {1/T:)5{r'^) = S{p)S{q). 

h—*o 

(4.36) 

Such a result is more or less to be expected since (4.30) 
shows that, for fixed n, En ^ as ?i ^ 0. 

A more interesting tack is to increase n as ?i — > so 
that En retains a constant value. One might imagine 
that then some classical features should emerge as ?i ^ 0. 
Without loss of generality, the energy may be taken to 
have the value 1/2, 

E = l/2. (4.37) 

So doing simply sets the scale in phase space. From (4.30) 
and (4.37) it follows that h is then given in terms of n by 
the relation 

?i = l/(2n-|- 1). (4.38) 

Correspondingly, the Wigner function takes the form 

W{n;r) = (-l)"(l/7r)(2n + 1) exp[-(2n + 1)7-^] 

xL^°\{4n + 2y]. (4.39) 

Classically the Liouville density is concentrated on the 
unit circle r = 1 when E = 1/2. Figures 1 through 3 
display W{n;r) as a function of r for n = 10, 20, and 
40, respectively. Evidently the ?i — > (n — > oo) limit 
in this case is far from simple, and requires some expla- 
nation. However, two features are immediately evident. 
First, the Wigner function does decay rapidly to zero for 
r > 1; but, contrary to classical expectations based on 
the Liouville density, it does not decay for r < 1 . Instead 
it oscillates rapidly and even increases as r decreases. In- 
deed, since Ln ^ (0) = 1, at the origin the Wigner function 
has the value 

M^(n;0) = (-l)"(l/^)(2n + l), (4.40) 
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which approaches ±00 as n — > 00. Second, the number of 
oscillations in the interval [0, 1] increases linearly with n. 
Consequently, derivatives such as (d/dp) produce factors 
of n. Equivalently, derivatives produce factors of 
Therefore, as before, what appear to be quantum correc- 
tion terms in (4.24) actually diverge as ?i like {l/K) 
and higher powers of {l/K). 
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FIG. 1 The Wigner function for n = 10. 



W{20;r} 

ft 

ft ft ft A A A A 






p V V V 

v 




r 



0.2 0.4 0.6 0.8 1 
FIG. 2 The Wigner function for n = 20. 
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FIG. 3 The Wigner function for n = 40. 



What sense can be made of the highly oscillatory be- 
havior of W{n\ r) for r < 1? One answer is again to 
view the Wigner function as a distribution. The den- 



sity operator corresponding to an ensemble, all of whose 
members are in the same energy eigenstate is given 
by the relation 

p^\n){n\. (4.41) 

Apply (3.11) to the operator {HY for £ = and 1. Mak- 
ing use of (3.3), (3.4), (4.25), (4.30), and (4.41) gives the 
relation 

/■OO 

2iT / drrW{n-r){r^ /2f ^tr[{HY\n){n\] 
Jo 

= [(n + l/2)f,]\ £ = and 1. (4.42) 
Now put (4.38) into (4.42) to find the result 

F{n, e) = 2TT drr^^+^W{n; r) = 1, l^Q and 1. 
^0 

(4.43) 

It is tempting to conjecture that (4.43) holds for all posi- 
tive integer values of ^. However, some machinery beyond 
that already presented is required for such a calculation 
because if A and B are operators with Weyl correspon- 
dences A and B, then it is not generally the case that 
the Weyl correspondence of the product AB is the sim- 
ple product AB (even if A and B commute, and even 
if they are equal). Often there are fi dependent correc- 
tions. They arise because there is an ordering problem 
for mixed operator expressions of the form (3™P". The 
direct and inverse Weyl transforms (3.1) and (3.2), which 
were used to define the Wigner function, imply a certain 
symmetric ordering procedure |^ . 

By a suitable change of variable, (4.43) can be written 
in the form 

F{n,l) = (-l)"[2(4n + 2)^]-iG(n,£) (4.44) 

where 

/•oo 

G{n,£)^ dxexp{^x/2)x'^L^^\x). (4.45) 
Jo 

Like all classical polynomials, the Laguerre polynomials 
obey a 3-term recursion relation. In their case the recur- 
sion relation reads |^ 

= -nL^°U^) + {2n + 1)4°) (x) ~ {n + l)4j,(a:). 

(4.46) 

It follows from (4.46) that the G(n, £) obey the recursion 
relation 

G{n,£+1) 

= -nG{n --l,£) + {2n+ l)G(n, £) - {n + l)G{n +!,£). 

(4.47) 
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Correspondingly, by (4.44), the F{n,£) obey the recur- 
sion relation 

F{n,e+1) 

= { [n(4n - 2)7(4n + 2)^+^] F{n -!,£) + {l/2)F{n, £) 
+ [{n + l)(4n + 6)V(4n + 2)^+^] F{n +!,£)}. (4.48) 

Finally, repeated application of (4.48) yields for the first 
few F{n, £) the results 

f(n,0) = l, 
F(n,l) = l, 

F(n,2) = l + [l/(2n + l)2], 
F(n,3) = l + [5/(2n + l)2], 

F(n,4) = 1 + [14/(2n+ 1)2] + [9/(2n + l)'^)]. (4.49) 

In view of (4.38), the terms in (4.49) involving n are, in 
effect, ft dependent corrections. 

Evidently (4.43) does not hold for £ larger than 1. 
However (4.49) shows that, for small £ and n — > oo 
{h ^ 0), there is the limit 

F{oo, £)= \im F{n,£) = 1. (4.50) 

n — >oo 

Moreover, for large n (?i — » 0) the recursion relation 
(4.48) takes the form 

F(oo,^+l) 

= (l/4)F(oo, £) + (l/2)F(oo, £) + (l/4)F(oo, £) 

= F{oo,£). (4.51) 

Therefore (4.50) holds for all £. 

Suppose the integral on the left side of (3.11) is to be 
calculated for any A{q,p) that is polynomial in the vari- 
ables q,p, and suppose that W{q,p) depends only on the 
variabk^ r. Them thc^ integral may be factored into an- 
gular and radial parts, and a moment's thought reveals 
that each angular integration vanishes unless the accom- 
paning power of r is even. Consequently, only integrals 
of the form (4.43) occur in any such calculation. Polyno- 
mials in q,p form a dense set of functions. Therefore, in 
view of (4.43) and (4.50), W{n, r) has (in the distribution 
theoretic sense) the n ^ oo (?i ^ 0) limit 

lim W^(n;r) ~ (l/27r)(5(r- 1). (4.52) 

n — ^oo 

Based on (4.52), one might conjecture that, in the limit 
n ^ 00, most of the oscillations of W(n; r), all of which 
occiir for r < 1, should nearly average to zero when 
integrated over any smooth test function. To examine 
this hypothesis. Figures 4 through 6 display the product 
2TTrW{n; r) as a function of r for n = 10, 20, and 40, re- 
spectively. Evidently there is a positive peak near r = 1 
preceded by oscillations to the left. The area under this 
peak is found to be 1.2638, 1.2679, and 1.2704 for n = 10, 
20, and 40, respectively. If this peak along with the full 
oscillation preceding it is considered, the areas are found 



to be 1.1387, 1.1460, and 1.1503 for n = 10, 20, and 40, 
respectively. Evidently the first few oscillations do nearly 
"saturate" the relation 

/•OO 

F(n,0)=/ dr2TTrW{n,r) = 1. (4.53) 
Jo 

However it is also evident that, even in the n ^ oo limit, 
it is necessary to retain all oscillations to fulfill (4.53) 
exactly. 
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FIG. 4 The integrand in (4.53) for n = 10. 
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FIG. 5 The integrand in (4.53) for n = 20. 
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FIG. 6 The integrand in (4.53) for n = 40. 
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V. DISCUSSION 

It has been shown that in general the Wigner function 
does not obey (1.11) for nonhnear symplectic transfor- 
mations, but rather there are "quantum corrections" as 
iUustrated by (4.24). This fact is no surprise to those 
well acquainted with Wigner function literature Let 
H he a, Hamiltonian of the form T + F as in (3.4), and 
assume that T is quadratic in p. Then it is known that 
the Wigner function obeys the equation of motion |^ 

dW/dt = -[W,H] 

-{h'^ /24)[{d/dqfH][{d/dpfW] + ■■■. (5.1) 

Retaining only the first (Poisson-bracket) term on the 
right side of (5.1) is equivalent to asserting (1.11). Note 
also that if H contains only linear and quadratic terms, 
which is equivalent to considering only translation and 
linear symplectic maps, then the correction terms in (5.1) 
vanish, and correspondingly (1.11) is exact. Suppose H 
is of the form 

H^-{a/3)q^ (5.2) 

The first correction term on the right side of (5.1) then 
gives the result 

- (h'^ / 2A) {d / dqf H = an'^/ 12, (5.3) 

in agreement with the first correction term in (4.24). 

It is also well documented in the literature that Wigner 
functions are often highly oscillatory, and that the ?i ^ 
limit is subtle . Correspondingly, it is well recognized 
that the correction terms in (5.1), just as the correction 
terms illustrated in (4.24), must be treated with care 
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